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ABSTRACT 

The published Mount Wilson Doppler-shift measurements of the solar velocity field 
taken in 1967-1982 are revisited with a more accurate model, which includes two terms 
representing the meridional flow and three terms corresponding to the convective limb 
shift. Integration of the recomputed data over the visible hemisphere reveals significant 
variability of the net radial velocity at characteristic time scales of 0.1-10 years, with 
a standard deviation of 1.4 m s^^. This result is supported by independent published 
observations. The implications for exoplanet detection include reduced sensitivity of the 
Doppler method to Earth-like planets in the habitable zone, and an elevated probability 
of false detections at periods of a few to several years. 

Subject headings: Sun: activity — planets and satellites: detection — techniques: radial 
velocities — stars: individual (55 Cnc) 



1. Introduction 



The spectacular success of the exoplanet search program, resulting in the discovery of over 300 
planets and planetary system to date, has been achieved mostly through the Doppl er-shift tech- 



nique . Indeed, spectroscop ic detections constitute the majority of known systems (jButler et al 



20061 : lUdrv k. SantosI 120071 ) . with the bias toward short-period, massive planets probably due to 
the selection effects of this method. Undoubtedl y, there is much room for further progress with 
the Doppler technique (jEggenberger &: Udryll2010l ) , with the accuracy of spectroscopic instruments 
steadily improving and now reaching ~ 1 m s~^ and the sensitivity of telescopes extending toward 
fainter stars. Combining precision photometric observations with radial velocity measurements 
provides a range of important physical characteristics of transiting stars invaluable for our under- 
standing the physics and the origin of exoplanets. The strategic goal of detecting rocky, habitable 
planets outside the Solar system now seems to be coming within reach. As the instrumental pre- 
cision steadily improved, a growing attention has been paid to the intrinsic perturbations in the 
observable parameters used in exoplanet detection. Stochastic, uncorrelated physical perturbations 
increase the level of noise, making it difficult to achieve the threshold signal-to-noise ratio, while 
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possible cyclic processes in the host s tars can mimic exoplan et signatures. In particular, the ro- 
tating pattern of photos pheric spots (ISaar &: Donahud 119971 ) and irregularities in the convective 
structure on the surface ( Meunier et al. 20ld ) can lead to an intrinsic scatter of radial velocities of 
up to a few m s~^ for solar type G dwarfs. In the younger Hyades, which are more magnetically 
active and rot ate faster than t he Su n, this effect is magnified to ~ 16 m s~^ in standard deviation. 



as observed by lPaulson et al.l (120041 ). As the level of activity is not constant in solar- type stars, the 
intrinsic radial velocity scatter is correlated with the magnetic cycle, opening possibilities of more 
sophisticated spectroscopic analysis in order to mitigate these difficulties. For example, the index 
of chromospheric activity is correlated with the area of star spots (a nd hence, with the observed 
scatter) and can serve as an indicator of magnetically induced cycles (jSantos et al.l l2000l ) . Careful 
selection of target stars can further improve the prospect of detection of smaller planets with the 
Doppler technique. Recent simultaneous measurements of chromospheric activity and radial veloc- 
ity imply that K-type dwar fs may be significan tly less variable than the Sun, and the intrinsic RV 



jitter is less than 1 m s ^ ( Santos et al.ll2010l ). Some stars ol der than 6 G yr and evolved off the 



main sequence have sharply reduced levels of magnetic activity (jWrightll2004l ) compared to the Sun. 
In this paper, we consider another important source of intrinsic radial velocity variation, related 
to the physical motion of the surface layers of stars, which remained largely outside the scope of 
previous papers. 

The surface layers of the Sun, where the spectroscopic lines are formed, are known to be 
involved in a complex pattern of radial and tangential motion. Furthermore, it is now an observa- 
tional fact that this velocity field is not static. In this paper, the series of Doppler measurements 
taken at the Mount Wilson Observatory is revisited, and the published fitting model coefficients 
are transformed to a different model of differential rotation, convective blueshift and meridional 
flow, which may more adequately represent the reality (§[2])- The transformed model is integrated 
to produce the net radial velocity of the Sun as a star (§ [3]). The results are compared to other 
data on the variability of the main components of the velocity field in § [U and a good agreement 
is found. The impact of the intrinsic variation of solar RV on the detectability of exoplanets, es- 
pecially within the habitable zone, are investigated in § [5] by means of ^-nd spectral density 
periodograms, as well as by a planet detection experiment, which results in two bogus planets. The 
relative importance of surface flows compared with other sources of intrinsic RV perturbations are 
discussed in § [H 



2. Mount Wilson data 



An extensive set of Doppler velocity measur ements taken at Mount Wilson during a period of 14 
years, from 1967 through 1982, was published bv lHoward et al.l (jl982l ). The full-disk magnetograms 
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m Fe I line 5250.2 A were fitted with the models] 

V = {A + B sin^ (j) + C sin^ (p) sin A cos (f) cos bo + 

{D + E sin^ </> + F sin'' (p) sin^ A cos (j) cos bo + G (1) 

with (j) the solar latitude, A the longitude (heliocentric angle from the central meridian), bo the solar 
latitude of disk center, G an arbitrary zero-point velocity, and A, B, C, D, E and F free fitting 
parameters. The first group of terms with coefficients A, B, and C represents differential rotation, 
with A denoting the equatorial velocity. Since the terms of differential rotation are odd functions of 
A, their net contribution to the integral radial velocity is zero, and we are not concerned with them 
in this paper. The second group of terms with D, E and F can be interpreted as the antisymmetric 
component of the differential rotation. These terms are even in A and result in a net contribution to 
the radial velocity. However, because the observations are differential and the term G is arbitrary, 
the total offset of velocity is indeterminate, and only the variation of these terms in time can be 
considered. 

The integral radial velocity as measured by a distant observer is 

ZL — 

- \ j 2 j 2 22 

V = — / d\ I {ao + Q-i cos (j) cos A + Q2 cos (j) cos A) V cos (j) cos A cos (jjdcj). (2) 
27r J__ir 

2 2 

In this formula, qq, ai and Q2 are the coefficients of solar limb darkening. Adopting = 0.30, 
ai = 0.93 and 03 = —0.23, the resulting net radial velocity is 

V = 0.1005 D + 0.01841 E + 0.007417F. (3) 



Using the tabulated data in (jHoward et al.lll982l ). we compute a radial velocity curve, which shows 
significant variations of 4.95 m s~^ in standard deviation (STD). However, this result is likely 
overestimated, for the following reasons. 



As was pointed out by (jSnodgrass &: Ulrichlll99d : ISnodgrass &: Howardlll985l ). the terms of dif- 



ferential rotation A, B s\v? (j) and C sin'* (p are not orthogonal on a half-sphere, resulting in consider- 
able correlations between the coefficients determined from regularly sampled data. They suggested 
to use Gegenbauer polynomials in sm(j) to orthogonalize the terms. The first thus orthogonalized 
coefficient corresponding to the principal component of solar rotation, \s A = A + B/5 + 3C/35. 
Note that the coefficients in this relation are closely in the same proportion (1 : 0.2 : 0.085) as the 
coefficients in Eg. IIOI (0.1005 : 0.01841 : 0.007417). This coincidence may be interpreted as a real 
net variation of the integral velocity (which can not be captured by the odd A, B and C terms) 
translating into the non-orthogonal D, E and F terms according to their covariances, since all the 
differential measurements are brought to a single zero-point. This hypothesis is supported by the 
actual negative correlation of the three terms in Eq. [T0\ whose variances add up to a larger number 
than 4.952. 



^The instrumental drift term has been omitted 
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Fig. 1. — The velocity field of solar meridional flow. The net motion is directed from the equator 
toward the poles, vanishing at the equator and the poles. 
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Thus, a real variation of the surface velocity field was probably exaggerated in the model 
([1]) because the DEF terms are internally correlated. Besides, these terms poorly describe the 
real velocity field. Indeed, the functional form of DEF terms corresponds to a pattern of zonal 
stretches on the surface, which remain fixed with respect to the central meridian despite the general 
rotation of the Sun. Clearly, any longitudinal perturbations of rotation velocity can not remain 
static with respect to the observer on time scales longer than one rotation period. The correct 
physical explanation to the symmet ric pattern of Doppler velocities observed at Mount Wilson was 
given by iLaBonte &: Howard! (jl982l ) . The mysteriously looking zonal pattern of radial velocity is 
caused by two separate physical phenomena on the Sun: the axisymmetric convective blueshift 
(often called the limbshift) and the meridional surface flow symmetric around the equator to first 
approximation . 



The model of meridional flow considered by iLaBonte fc Howard! (119821 ) assumed a constant 
surface velocity directe d toward the poles. We replace it with a more adequate model suggested by 
Hathawav et al.l (|l99d ): 



Vm = /i2-F2 (sin 0) + /U4P|(sin 



(4) 



where are associated Legendre polynomials of degree I and order 1. The free fitting parameters 
^2 and /i4 are to be determined from observations. The Legendre polynomials are orthogonal on the 
range of sin (j), which drastically simplifies the subsequent analy sis. We also repl ace a polynomial 
expansion for limbshift with the more sophisticated model from (jHathawavl Il996l ) : 



Vb 



k=l 



ftP;*(l-COS p), 



(5) 



where are shifted Legendre polynomials orthogonal on [0,1], p is the central angle (cos p 
cos (j) cos A), and /3j are free fitting coefficients. 

Thus, the new model, which more faithfully represents the observed velocity field is 



V = {A + B sin^ 
{P2P2 (sin 4 



> + C sin^ (f)) sin A cos <f> cos 60 + 



■ P4-P4 (sin (^))(cos A sin (f) cos 60 — cos (j) sin ho) + Vb + G. (6) 



We carried over the ABC terms of differential rotation (rather than using Gegenbauer polynomials) 
to simplify subsequent transformations. 

Our task is now to recompute the observed parameters of the velocity field using the published 
ABCDEF parameters. The fitting procedure is a linear least-squares problem and the transfor- 
mation problem can be expressed in matrix notation 



[A i = [A T] X + e. 



(7) 



where A is the matrix of differential rotation (the ABC terms), $ is the matrix of the symmetric 
stretch (the DEF terms), T is the matrix of the new meridional flow and limb shift terms, x is 
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the vector of free parameters of the old model, x is the vector of free parameters of the new model, 
and e is random noise. It can be shown that in a least-squares solution, 



X = TT [A X, 



(8) 



where = (x"'"T)~^T"'' is the pseudoinverse of T. The matrix and inner vector products here 
are calculated through integration of the element products over the visible disk with the cos b 
weight. The inner products of the new model terms and the ABC terms are all zero (in other 
terms, these functions are orthogonal). Therefore, the ABC coefficients remain unchanged in this 
transformation, and the new parameters depend only on the DEF coefficients. 

After some toil with integration and matrix inversion, one obtains 



^^2 




-0.766 


-0.252 


-0.119 


/"4 




0.482 


0.094 


0.018 


Pi 




0.345 


0.108 


0.053 


P2 




0.012 


0.030 


0.026 


/33 




-0.106 


-0.030 


-0.010 



D 
E 
F 



(9) 



3. The solar RV variation 

Analogous to Eq. [2l the new model of the solar velocity field can be integrated over the visible 
hemisphere, to yield: 

V = -0.157/i2 +0.025 /i4 - 0.295 /3i +0.044/32 +0.004/33. (10) 

This result is valid only for the Sun seen edge-on (inclination 90°), although the limb shift is 
probably independent of inclination. The five terms of the new model are derived from only three 
terms of the old model, hence they are strongly correlated because of the projection i n the function 



space . The reconstructed radial velocity curve is shown in Fig. [2j The original data in ([Howard et al 



19821 1 are sampled by Carrington rotations, therefore, each data point represents the mean over 
27.2753 days. Monthly variations of a few m s~^ are common, but so are longer-term variations on 
time scales of years. A running median shown with the red line indicates that the long-term radial 
velocity of the Sun grew by ~ 7 m s~^ between 1976 and 1981. The reconstructed RV curve has a 
standard deviation of 1.44 m s~^ . The long-term behavior of solar RV is strongly correlated with 
the magnetic cycle (cf. the monthly sunspot numbers in the lower panel of Fig. [2]). Radial velocity 
variations on time scales of several months to 1 year have a strong impact on the detectability of 
habitable planets around solar-type stars, as discussed in § [5l We have to carefully investigate if 
these results are consistent with other independent data on the Sun. 
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Fig. 2. — The integrated radial velocity of the Sun seen equator-on reconstructed from the Mount 
Wilson Doppler measurements (upper panel). A running median of 7 individual data points (one 
per Carrington rotation) is drawn with the red line. The zero-point is arbitrary. In the lower panel, 
the sunspot number during the same period is shown for reference. 
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4. Comparison with other data 



Noting that the low-order terms of meridional flow (fj,2) and limb shift propagate with 
the largest weights into the integral velocity in Eq. [TOl we have to consider what is known about 
the variability of these physical effects. The most direct comp arison comes from the Doppler 
observations with the Global Oscillation Network Group (GONG). iHathaway et alJ (119961 ) present 
the results for the main component of the meridional flow over some 130 days in 1995 (their Fig. 1, 
top panel), corresponding to P2 (Legendre polynomial of degree 2 and order 1). These coefficients 
are denoted ^2 in Eq. [TUJ The flow velocity varies between and 50 m s~^, which for an edge-on 
Sun translates to a peak-to-peak amplitude of 8 m s~^. Quasi-periodic cycles of 10-20 days are 
evident in the GONG data, as well as longer-term variations on a few months, which are expected 
to change the integral velocity by roughly 3-4 m s~^. The amplitudes of the higher-degree terms 
P4 and Pq are found to be much smaller (—1.6 it 0.5 and —1.3 it 0.1 m s~^), so their contribution 
to the net radial velocity is negligible in the context of this paper. 

On lo n ger ti me scales, the main I = 2 component of meridional circulation was investigated by 
Hathaway! (|l996l ). also based on the GONG Doppler measurements. The flow velocity grew from 
~ 20 m s~^ in the summer of 1992 to almost 100 m by the end of 1993, then dropping to 
small negative values in the summer of 1994, indicating a brief episode of reversal. After that, it 
gradually returned to ~ 20 m s~^ by the mid 1995. From Eq. 1101 the overall radial velocity should 
have changed by 16 m s^^ in the first half of 1994. Unfortunately, our data in Fig. [2] cover an 
earlier time interval, and can not be directly compared with the GONG results. But the scale of 
variability in our reconstruction is absolutely consistent with the newer, and more accurate data. 



In the same paper, iHathawayl (|l996l ) presents a reconstruction of the history of the three terms 
of convective blueshift, denoted (32 and in this paper. Only the first term appears to be a 
significant contributor to the net velocity, with a quasi-sinusoidal variation of ~ 30 m s~^ on a 
period of 1.5 yr and a short-term scatter of similar magnitude. The corresponding range of radial 
velocity is 0.295 • 30/^/3 = 5 m s~^. Little correlation was noted between the meridional flow and 
the hmb shift on the GONG data. 

The meridional flow traced by the motion of magnetic features in the photosphere from a 
26-year set of Mt. Wilson magnetograms appears to be somewhat smaller in amplitude and more 
complex in the equatorial zone, with strong equator- ward features (ISnodgrass &: Dailevlll996l ). The 



maximum rate is ~ 13 m s~^, and the pattern is strongly correlated with the solar cycle. A 
det ailed picture of the sola r cycle-related variation in the structure of the meridional flow is shown 
in (jUlrich &: BovdenI |2005| . their Fig. 2). The poleward components of the meridional circulation 
are the stronges t when the solar a ctivit y is at its maximum, in agreement with the reconstructed 
data in Fig. [2l Gizon Sz Rempel (j2008l ) used MDI full-disk Doppler images for the period 1996- 
2002 and measured the advection of the supergranulation pattern, which is probably consistent 
with the tracing of magnetic features. They detect only poleward motions reaching 13 m s~^ 
in the antisymmetric part of the field averaged in 1 yr intervals. The peak-to-peak variation of 
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the antisymmetric component is 7 m s ^. This is again in agreement with the magnitude of RV 
variation derived in this paper. 

With the abundance of high-quahty data, the analysis procedures are very complex, which 
accounts for some lingering disagreement between different published result s on the so l ar vel ocity 



field. One of the most rigorous and comprehensive models was presented by iHathawayi (|l987l ). in- 
cluding differential rotation, torsional streams, meridional circulation, the limb shift, supergranules 
and giant cells. The surface velocity field is represented with vector spherical harmonics, which 
are orthogonal on a unit sphere. Although this much desired orthogonality is compromised by 
the availability of only one hemisphere, inclination and projection effects (as far as Doppler data 
are concerned), the various scales of poloidal and toroidal motions are naturally separated in the 
model through vector spherical harmonics of different degrees and orders. For example, the three 
main components of differential rotation are represented by the toroidal harmonics of zero order 
(solid body rotation), Tg and Tg. The limb shift terms are correlated (non-orthogo nal) to the 



merid ional flow terms to such extent, that their separation requires a full-cycle iteration (jHathaway 



19921 ) ■ especially if the aim is to reconstruct the absolute velocity field. This seems to be the only 
way to infer the main term of limb shift (the value at the center of the disk), which is —540 m s~^ 
according to Hathaway. Thus, the variable component of the limb shift is only w 1% of the constant 
component. 

Similarly, the large-scale meridional flows are dwarfed by t he supergranulation pattern, which 



has typical velocities of 300 to 400 m s ^ ( Hathaway et al.ll200d ). The high-degree supergranulation 



is of little interest in the context of this paper, because it only contributes to the high-frequency 
variation of the integral radial velocity, whereas we are interested in what happens with solar radial 
velocity on the time scales of months to years. But the lower-degree supercells are of interest, 
because they rotate with the Sun, indicating physical flows in the photosphere. A rotating spherical 
harmonic term of order m will produce a periodic signal in the integral velocity with a period of 
PQ,rot/'m- The power spectrum of velocity fleld is peaked at degrees I ~ 120, which correspond to 
the spatial scale of supergranules. Therefore, rotating supergranules will cause a net variation in 
radial velocity of a few m s~^ on time scales of a few hours. 



5. The impact on exoplanet detection 

Exoplanets are found in Doppler measurement data through the periodic perturbations caused 
by the orbital motion of the host star around the system barycenter. To a good approximation, 
the radial velocity of the star is the sum of a constant term (the velocity of the barycenter) and 
a number of Kepler motions. For low-eccentricity orbits, such as those in the Solar System, the 
Kepler velocity signal is almost perfectly sinusoidal. Detecting one or several sinusoidal variations in 
high-cadence Doppler measurements superimposed with a constant offset is fairly straightforward, 
if the observational noise is random, uncorrelated, and sufficiently small compared to the expected 
signal. 
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An additional source of noise or perturbation from a real physical process on the star makes 
this detection significantly more difficult, or in some cases, ambiguous. To quantify the impact of 
the 'surface flow noise' on the detectability of planets, a generic planet detection algorithm was 
applied to the data derived in §[3l No additional noise was added to the individual RV points. The 
procedure starts with a periodogram analysis and a subsequent estimation of the confidence of the 
most prominent sinusoidal variations. Fig. [3] shows the resulting periodogram of the reconstructed 
solar radial velocities. Instead of the traditional spectral power, the residual is plotted for a 
high-cadence grid of periods. The statistics is preferred, because the confidence level of the 
detected minima is rigorously computed via the cumulative probability function of the F-statistics 
with {nobs — 3, 2} degrees of free dom. This confid ence F-test is analogous to the estimation of 



FAP (false alarm probability) bv ICummina (|2004l ). with the number of independent frequencies 



corresponding to the frequency search interval. Each dip in this curve which corresponds to a value 
of F surpassing a given confidence threshold, can be considered a planet detection. Setting the 
confidence threshold at 99%, this periodogram analysis results in the detection of two bogus planets 
with orbital periods Pi = 9.35 yr (mass 26M0) and P2 = 6.35 yr (mass I5M0). The former false 
detection is apparently triggered by the long-term variation of the meridional flow with the solar 
cycle, and the latter may also be related to some poorly understood cyclicity of magnetic activity, 
because a sinusoidal variation of similar frequency emerges in the Total Solar Irradiance (TSI) data. 

This numerical experiment indicates that the magnetic cycle-related variations in Sun-like stars 
can in some cases be confused with long-period planets. It is therefore of interest to compare the 
reconstructed solar periodogram with accurate RV data for a planet host. Fig. [H sho ws a peri 



odog r am for the star 55 Cnc with high-resolution spectroscopic measurements from (jMarcy et al 



2002; iFischer et al.l 120081 ). The total span of these measurements is 18.3 yr, and the formal sin- 
gle measurement error progressed from ~ 9 m s~^ during the first 6 years (Lick observations) to 
between 3 and 5 m s~^ over the later period (Keck observations). The star is a K0/G8 dwarf, 
slowly rotating and chromospherically inactive (logi?^j^ = —4.84). It is a close, albeit soin ewhat 



smaller, solar analog. Five planets have been reported in the literature ([Fischer et al.l 120081 ) . with 
orbital periods from 2.8 d to 14 yr. Fig. H] depicts a portion of the x^-periodogram computed 
from the weighted RV data for the same range of periods as Fig. [3l Our confidence estimation is 
consistent with the FAP thresholds computed in the cited papers, where a different technique was 
used (Monte-Carlo trials on randomly scrambled data). Except for a number of pronounced high- 
frequency dips, some of which correspond to the detected planets, the structure of the periodogram 
is quite similar to that for the Sun. In particular, the lowest dip in this part of the periodogram 
centered on 14 yr at the long-period end, corresponding to the planet 55 Cnc d, is reminiscent of 
the 9.35-yr feature for the Sun, which is quite likely caused by the solar cycle. This similarity calls 
for caution in the interpretation of long-period signals from stars like 55 Cnc. 

Extrapolating the observational results presented in § [2] to other stars requires caution. The 
Mount Wilson data are based on a single Fe line, whereas an RV curve for exoplanet detection is 
derived from a large number of weak lines across the visible spectrum. There are good reasons to 
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Period, years 



Fig. 3. — x^'Psriodogram of the reconstructed solar radial velocity curve in Fig. [2j The dashed 
horizontal line indicates the threshold corresponding to a false alarm probability of 1% for the 
first planet detection. 
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assume, however, that the dominating //-terms of the meridional flow in Eq. [6] are not wavelength- 
dependent. The existing models of motions in the Sun (still somewhat u ncertain) i mply large-scale 
and deep-rooted dynamical structures of global character. For example, iHowardI (|l987l ) presented 
a model of only several gigantic longitudinal supercells (or rolls), which may evolve or drift with 
the magnetic cycle. In this interpretation, the large-scale meridional flows are quite deep. The 
case with convective limbshift, which is wavelength-dependent, is perhaps more complicated. The 
blueshift is not the same in different lines, but it seems plausible that the long-term changes due to 
the 11-yr and other known magnetic cycles are caused by a slower or faster convective motion on the 
integrated disk, thus resulting in a uniform shift of line bisectors across the sp ectrum. Systematic 
changes in line shapes can be expected due to magnetic cycles as well. Recently, ISantos et al.l (j2010l ) 
investigated a sample of eight late G or early K dwarfs of low magnetic activity and compared 
precision RV measurements with simultaneous estimates of chromospheric activity. They could not 
find any clear correlation and concluded that any RV variations induced by magnetic cycles on 
these stars are unlikely to be much greater than 1 m s~^ . Therefore, it is possible that many K 
dwarfs are less subject to the adverse effects of variable surface flows than more solar- like stars. 




Period, years 



Fig. 4. — x^'Psriodogram of the RV measurements of 55 Cnc, which were used to detect the 
5 exoplanets reported in the literature. The dashed horizontal line indicates the threshold 
corresponding to a false alarm probability of 1% for the first planet detection. 



Most of the currently known o r suspected exopla i iets have masses of t he Solar System giants 
(or larger), and short-period orbits (jButler et al.ll2006l : lUdry fc Santosll2007l ). A considerable effort 
is under way to improve the sensitivity of the Doppler instruments and data calibration to the 
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levels sufficient to detect Earth-like, rocky planets within the habitable zones of main sequence 
stars, where water can be present in the liquid state. The required single measurement error due 
to photon noise and instrument imperfection should be < 0.1 m to achieve this goal. One 
should consider the possibility that the surface flow and convective blueshift perturbations are in 
fact much larger than 0.1 m s~^ in the range of habitable orbits too. Using the reconstructed solar 
RV curve in Fig. [21 I computed the spectral power density (i.e., the amplitude of the best-fitting 
sinusoid as a function of frequency) for a range of periods 0.5 through 1.5 yr, roughly corresponding 
to the solar habitable zone. The power density turned out to be quite flat at ~ 1 m s~^. This 
result is of modest value because it only sets the upper bound for the actual power spectrum of 
RV variations of the Sun, since the contribution of observational error in the Mount Wilson data 
is not known. Ultra-precise, long-term RV observations of stable solar-type stars will be required 
to confirm that the intrinsic RV variations come up to such values. We note that a 3a detection of 
Earth in a periodogram requires the perturbation power density to be less than 0.03 m s~^. 



6. Discussion 



Variable surface flows are only one of the known physical processes on the Sun that can change 
the integrated radial velocity. They appear to be long-term in nature, and therefore, may be the 
main obstacle to detecting Earth-like habitable planets with ultra-precise Doppler measurements. 
The main pulsation modes of solar-type stars (p-modes) have periods of 10-20 min, and can be 
successfully suppressed in observations by taking longer exposures. The supergranulation pattern, 
discussed in §|31 should produce variations of up to a few m s^^ but their time scale is several hours. 
Rotation of stars and the non-uniform distribution of surface brightness due to photospheric spots 
and plages is probably the main source of RV perturbations on the time scales up to 50 d. The 
impact of spots has been investigated in numerous papers, and the lat est estimates suggest relatively 
small, but non-negligible , dispe rsion for the Sun of ctrv — 0-4 m s~^ ( Makarov et al.ll2009l ). As was 
shown by lPaulson et al.l (j200J) for the Hyades, the photometric and the RV jitter of more active 
stars than the Sun are strongly correlated. Extrapolating the empirical relation found in that paper, 
one would expect a jitter of cjrv = 1-7 m s~^ for the Sun. The large discrepancy between these 
estimates is probably related to the different morphology of photospheric inhomogeneities in the 
Hyades stars. Dark spots are the dominating magnetic features on active and young stars, whereas 
the contributio n of bright plages on o lder, solar-type stars tends to match the impact of spots, or 
even to prevail (jLockwood et al.ll2007l ). Besides, active stars often have only one or two giant long- 
lived spots on the surface, resulting in a strong rotational modulation of the light curve, whereas 
the surface of an active s un is usually rnarked with a few spot groups at a time, fairly uniformly 
distributed in longitude. iMeunier et al.l (120101 ) estimated the combined impact of sunspots and 
plages at ~ 1 m s~^ in RV amplitude, and concluded that it should be conflned to periods less 
than 100 d. Our study suggests that on the time scales 0.6-1.4 yr, characteristic of orbits within 
the habitable zones, slowly evolving surface flows and the distribution of the convective blueshift 
is the major source of confusion and error. 
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Our conclusions are in agreement with the observations of solar radial velocity bv lMcMillan et al 



(|l993l ). These authors observed the solar light reflected from the Moon for 5 years nightly in vio- 
let absorption lines and determined an upper limit of 4 m for the overall dispersion of radial 
velocities. The periodogram of RV variations (their Fig. 4) shows four peaks rising above the 3cr 
threshold, all in the long-period part of the spectrum. Three of them may be of instrumental origin, 
but the unresolved peak at period longer than 3 yr appears to be genuine and is called "intriguing" 
in the paper. The two bogus planets found in the recons tructed RV data in thi s paper (§ [5]) have 
periods 6.35 and 9.35 yr. Therefore, it is possible t hat [McMillan et al.l (119931 1 already detected 



long-term variations in the RV of the Sun as a star. iMeunier et al.l (|20ld ) pointed out that the 
effects of convective blueshift, which is included in our RV model (§[2]), may be smaller in the deep 
violet lines than in the redder lines that are normally used for planet search. 

Our final note is that the model of meridional flow depicted in Fig. [T] is a simplifled representa- 
tion of the actual fine structure of the velocity fleld. By necessity, the model for the Mount Wilson 
data is coarse, where the contribution of convective blueshift and meridional flow can not be decou- 
pled. The former is especially problematic, because it is wavelength-dependent, and uncalibrate d 
changes in the spectrograph setup can bring about long-term systematic errors (cf. UlrichI [2OO1I ). 
The low-order polynomials used to flt the surface velocity field (§ [2]) only reveal the underlying, very 
much smoothed structure of the actual distribution of tangential and radial motion, characterized 
by smaller spatial scales and larger amplitudes . For examp le, the local flows are known to converge 
toward active regions, reaching 50-100 m s^^ ( Gizonll2004l ). and these local inflows may be respon- 
sible for the temporal variation of the overall velocity fleld on the time scales of supergranulation (8 
hours) and the solar cycle (11 years). Thus, the temporal behavior of the velocity fleld may be in 
part the average of many sto chastic components, and in part the evolution of the global magnetic 
field (jUlrich fc: Boydenll2005l ). both still poorly understood for the Sun. Therefore, it is doubtful 
that a good diagnostics can be devised for other stars to separate the physical motion of the surface 
from the signatures of small planets. 



The research described in this paper was carried out at the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under a contract with the National Aeronautics and Space Admin- 
istration. 
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